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1. Introduction 

It is well known that the time development of reaction diffusion systems in which 
two types of particles A and B ('particles' and 'antiparticles') annihilate irreversibly 
is dominated by density fluctuations [I], @, 0, f§, [5|]. If an equal number of A and B 
particles is randomly placed on the sites of a d- dimensional lattice the densities are 
asmptotically given by n(t) ~ t~ d ^ (for d < 4) whereas a naive approach based on 
mean-field rate equations predicts n(t) ~ t -1 . 

The slow density decay is due to the asymptotic segregation of particles into regions 
of purely A or B particles. The reaction becomes more efficient if the effects of 
segregation are reduced by an appropriate mixing mechanism like turbulence |J or 
random forces in low- viscosity liquids . In these cases the diffusive particle motion is 
replaced by a superdiffusive behaviour characterized by a mean-square displacement of 
the form (x(t) 2 ) ~ t 2 ^ z with z < 2. Very recently the annihilation reaction in a driven 
diffusive system with repulsive interaction between particles of the same type 0, has 
been studied. If the particles are subject to an external driving force in one direction the 
repulsion gives rise to superdiffusion with z = 3/2 (for d — 1), and the density decreases 
asymptotically as n(t) ~ t -1 / 3 . The relation between the long-time behaviour of n{t) and 
the value of the exponent z has been investigated for particles performing independent 
Levy walks |TD| . In this case z is defined through the scaling of the transition probability 



(which is a function of r~ z t) since the the mean-square displacement for Levy flights with 
z < 2 diverges. It has been shown that below the upper critical dimension d c (z) = 2z 
the long-time behaviour of the density is given by n{t) ~ t~ d ^ 2z \ 

The mixing caused by the inhomogeneous velocity field in a linear shear flow 
reduces the critical dimension to d c = 2, i.e. n(t) ~ t~ x in any physically 



accessible dimension [|ITJ. Other models for particle mixing have been considered in 
references |12, In this paper we investigate the A +B- annihilation reaction in a 
medium with a random velocity field. It is assumed that the velocity at every point 
r = (x, y) of a (^-dimensional system is parallel or antiparallel to the x-axis and depends 
only on the coordinates perpendicular to the flow. The velocity field is modelled by 
quenched Gaussian random variables with zero mean and the correlations 

KyMy')] = /*(y - y') (i) 

where the square brackets indicate the average over the realizations of the flow. 
Originally this model was introduced to describe the ground water transport in 



heterogeneous rocks |14] 



Below three dimensions a random walk in the presence of a velocity field of the 
form (TJ) displays superdiffusive behaviour in the x-direction [JT^, nj. For a particle 
starting at r = at time t — the mean square displacement in the x-direction averaged 
over the configurations of v(y) is given by [{x 2 )} ~ a 2 t^ d ^ 2 for d < 3 (with a generalized 
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diffusion constant a 2 ). One can use this result for a naive estimate of the densities at 
large times t. If we approximate the motion of the particles by independent anisotropic 
Levy flights with exponents zii = 4/ (5 — d) (in x-direction) and z± = 2 & straightforward 
generalization of the arguments given in |10| yields 

n A ,B(t) ~ [n^^U-^l^] 112 = nl /2 t-^ 8 (2) 

where Uq is the initial density of each particle type. However, it is not clear in how far 
the effects of the velocity field can be described by independent Levy walks since v (y) is 
a quenched random variable with an infinite correlation length in the x-direction. In the 
following sections we will use the field theoretic methods developed in [1I7L 0| to show 



that the power law (0) is correct for t — > oo and calculate the amplitude at first order 
in e = 3 — d. The result of this calculation reads 

n A , B (t) ~ A(e) (2vr 3 / 2 (87r J D)( d - 1 )/V) _1/2 yfarW* (3) 

where 

A( e ) = l + l(31n2-l) + ( e 2 ) (4) 

and n = n A (0) = tib(O) is the initial density of each particle type. In contrast to the 
case of annihilating Levy particles which can be treated in a mean field approximation 
with random initial conditions the correlated particle motion in a quenched velocity 
field requires the application of the renormalization group. 

In the next section a microscopic model for a reaction diffusion system in a random 
velocity field is defined and its mapping to a continuum field theory is discussed. In 
Section |] we briefly review the renormalized field theory for non-interacting particles in 
a random velocity field. The results are used in Section |] to derive an effective action 
which describes the long-time behaviour of the system. The asymptotic decay of the 
densities is calculated in section || and the results are discussed in Section ||. 



2. The model 



At every time t the microscopic state of the system is defined by the occupation numbers 
m(rj) of A-particles and n(r^) of B-particles at every site i (with position r») of a d- 
dimensional lattice. Denoting the probability for a given configuration {m, n} at time t 
by P({m, n}; t) we can describe the stochastic dynamics of the system by a set of linear 
differential equations (master equation) for the probabilities P({m,n};t). (We assume 
that the time evolution is Markoffian.) In order to apply field theoretic renormalization 
group methods we use a functional integral description which is equivalent to the master 
equation. Since a detailed derivation of this formalism is given in references [O, |1^. Jl 
we only give the main steps focussing on the modifications that are necessary to study 
the influence of the random velocity field. 



We first map the probabilities P({m,n};t) to a vector of the infinite 

dimensional Fock space spanned by the vectors \{m,n}) by writing |E| |2(J 

!*(*)>= E P({m,ny,t)\{m,n}). (5) 

{m,n} 

Since the master equation is linear and of first order in t it may be written in the form 

d t mt)) =imt)) (6) 

with an appropriate Liouville operator L. This operator can be expressed in terms of 
the annihilation operators aA(r«), flB( r i) and creation operators a\(ri), a B (rj) defined 
by 

a A (r i )|m(r i )) = m(Y i )\m(r i ) - 1) a\(T^\m{T$) = |m(rj) + 1) 
o B (ri)|n(ri)) =n(r i )|n(r i ) -1) ^(r^Mr*)) = {nfa) + 1). 

These operators leave the occupation numbers m(rj) and n(rj) of the sites j ^ i 
unchanged. 

The diffusion of non-interacting particles corresponds to the Liouvillean 
^ = t|E («A( r i) - OaO*)) (oA(rO - a A (Tj)) 

11 <i,j> 

+T5 E ( a B( r i) " « B ( r i)) (obW - aB(rj)) (8) 

n <i,j> 

where the sum extends over all pairs of nearest neighbour sites, h is the lattice spacing 
and D denotes the diffusion constant. The annihilation of particle-antiparticle pairs at 
the same site is described by 

k 

^reac = E i 1 " S A ( r i)«B ( r i) ) « a(* 'OM 1 \) (9) 

i 

where k/h d is the reaction rate. 

To investigate the effect of a random shear flow we introduce independent Gaussian 
random variables v(y) labelled by the coordinates y perpendicular to the direction of 
the flow. At each point rj = (xj, y«) the velocity field prefers particle jumps parallel 
or antiparallel to the x-direction depending on the sign of v (y^). If v(yi) is positive 
particles at the point (xi, y^) will jump with rate f(y;) jh in positive x-direction whereas 
for v(yi) < they will hop with rate —v(yi)/h in the opposite direction . Here we 
assume that v(y) have zero mean and the correlations 

Hy)v(y')\ = fh % y . (10) 

While deviations from the Gaussian distribution turn out to be irrelevant for the 
asymptotic scaling behaviour long range correlations of the velocity field [|2T| or a non- 



vanishing mean shear 1TTJ change the universality class. The motion of the particles in 



the velocity field is described be the Liouvillean 
1 
h 



iflow = T^2\v(ji) \ [a A (r» + sgn(>(y i))he x ) - a A (r;)] a A (rj 



+\l2\ v (yi) \ l d B( r i + sgn(v(yi))he x ) - ^(r»)] asfo) (11) 

where sgn(f ) = ±1 is the sign function and e x denotes the unit vector in positive 
^-direction. In order to perform the average over the realizations of v(y) it will be 
convenient to write Lfl ow as a sum, Lfl ow = + L^ + \ where 

= Yj l^l^i + he x) - a A ( r i ~ he x)] 3A(r») 



i 



+ 2^£ v (y<) l^B^i + he x) - a B ( r i - he x )]a B (ri) (12) 



is odd with respect to the velocity field and 
1 

2h 



\ v (yi) \ + he x) + a* B (ri - he x ) - 2c&(r i )] a B (r<) (13) 

i 

depends only on the modulus |f(y)|. The even part can be interpreted as a 
contribution to the diffusion in x-direction with a y-dependent diffusion constant. 
The formal solution of equation (j^j is given by 

|$(t))=exp(Lt)|$(0)> (14) 

with L = Ldiff + I/rcac + iflow In order to derive a functional integral representation for 
the dynamics one uses the Trotter formula 

exp(Lt) = lim (1 + Lt/n) n (15) 

n— >oo 

to rewrite the time evolution operator as a product of n factors linear in L. Inserting the 
identity operator in a coherent state representation between the individual factors jLT], 
|l~8l , U one obtains a path integral with action 

S[a, a; b, b] = h d J dt ^[a(r f , t)d t a(r h t) + 6(r 4 , t)d t b(r h t)] 

i 

- J dtL[a(t), a{t); b{t), b{t)\ - h d £ n (a A (r 4 , 0) + a B (r 4 , 0)) (16) 

i 

where a, a, 6 and 6 are c-number functions and L[a(t),a(t);b(t),b(t)] can be obtained 
from L via the replacements 

a A (ri) -> /i d a(r i; t) a A ( r i) 1 + o( r ^ *) ^ 
a B (ri) -> h d b(r u t) a^fa) -> 1 + &(r i5 t). 
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In equation (16) n denotes the initial density of each species. 

In order to calculate the average of correlation and response functions with respect 
to disorder we use the effective action S defined by 

exp(— S[d, a; b, b}) = exp(— S[a, a; b, b]) . (18) 

For each value of the coordinate y one has to calculate an integral of the form 

My)(2n ///^V^exp (^L) ' «p(-At;(y) - B\v(y)\) (19) 

where A and B are shorthand notations for the contributions to S coming from L^~' 
and L^ + \ respectively. The integral (|T9|) can be expressed in terms of error functions. 
We only need the approximation 



m = exp ( lif/h^A* - ^{j/h^f^B 



(20) 

since higher orders in A and B turn out to be irrelavant for the asymptotic scaling 
behaviour. The first term in the exponential function (proportional to A 2 ) gives rise to 
a new interaction which is non-local with respect to time while the second contribution 
(linear in B) modifies the diffusion constant in the direction parallel to the flow. 
Neglecting all irrelevant terms one arrives at 



S[a, a; b,b] = J dt J d d r [a(d t a - DA ± a - D\\d\a 



+ b(d t b - DAJ - D\\dp) + kab(a + b + tib)- n (a + b)S(t) 

~ 2 I ^ ^ (/ dt / dx ^ 9 W a + ^H 6 )) 2 ■ ( 21 ) 
Here the sum over lattice sites has been replaced by an integral over the continuous 
variable r = (r», r_|_) = (x, y). This continuum model is appropriate for the study of the 
long-time and large- distance behaviour of the system. 

Since the quantity ip = (a — b)/ y/2 is closely related to the conserved density 
difference of A-particles and B-particles it is convenient to introduce the fields Q 

i, = -L(a-b) ^=-L(5-6) <p=^=(a + b) = -L(o + 6). (22) 

After rescaling of time by the diffusion constant D the action becomes in terms of the 
new variables 

S$, & 4>, 4>] = J dt J d d r [i>(d t ij - A ± ij - Kd\i)) + 4>{d t cj) - A ± 4> - Kd\<P) 
+\ l 4>(<j ) 2 - ^ 2 ) + \ 2 (4> 2 - ^ 2 )(0 2 - ^) - V2n 4>5(t)] 



d d - l y (j dt J dx^d^ + (23) 
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with the coupling coefficients Aj = k/(y/2D), A 2 = k/(4D), f = f/D 2 and k = D\\/D. 

The action S includes two different types of fluctuations. The interaction 
proportional to A2 leads to anticorrelations neglected in the mean field rate equations. 
In the case of the single species annihilation A + A — ■> these anticorrelations are 
responsible for the slow density decay n(t) ~ t~ d l 2 in d < 2 [18|. The interaction 



proportional to / describes the effect of the quenched random velocity field and gives rise 
to the super diffusive motion of particles found in references fi"6| , jl5f . A straightforward 
dimensional analysis shows that A2 is relevant below d\ = 2 while the upper critical 
dimension of the disorder is df = 3. As already pointed out in [|J this does not mean 
that A2 may be neglected from the outset if d > d\. For no n- vanishing A2 coarse graining 
of the action S generates new interactions which are relevant below four dimensions. 
Since these interactions are located at the 'time surface' t = they modify the initial 
distribution of the fields ip and 0. An effective action appropriate for the analysis of 
fluctuation effects in dimensions d > d\ includes all relevant interactions which may be 
generated by the renormalization group for A2 > 0. 

We will use renormalization group improved perturbation theory to derive results 
from the effective action for non-zero /. Since this calculation is based on the results of 
reference JOj] for non-interacting particles (i.e. Ai = A 2 = 0) it is useful to review the 
renormalization group approach to this problem briefly in the next section. 



3. Non-interacting particles in a random velocity field 

For Ai = A2 = no = only Green functions of the form 

I M M+N 

G MtN {{r, 0) = ( II ^(r*, U)4>(ii, U) ■ J[ 0(r i; %)fe> *i) 

\i=l j=M+l 



(24) 



A=0 



with an equal number of insertions of ip (0) and if) (0) are non-zero. (Throughout this 
section the angular brackets indicate an average with weight exp(— S) for Ai = A2 = 
n = 0.) 

Since / is a relevant parameter below three dimensions the asymptotic scaling 
behaviour of the Green functions cannot be obtained from a naive perturbation 
expansion at finite order in /. We therefore use the renormalization group to improve 
the perturbation theory by a partial summation of the perturbation series. Due 
to ultraviolet-divergencies at the upper critical dimension df — 3 the individual 
contributions to the perturbation series have poles in e = 3 — d. In the minimal 
renormalization scheme these poles are absorbed into renormalizations of coupling 
coefficients and fields. 
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The propagator of the free (Gaussian) field theory is given by 

(V>(r, t)${0, 0)) = (0(r, t)4>(0, 0)) = / e^ r G(q, t) (25) 

with the Fourier transform 

G(q, t) = exp (-(<£ + «gf)t) . (26) 

Here B(t) is the Heaviside step function and the indices '||' and indicate the respective 
directions parallel and perpendicular to the flow. Due to the strong anisotropy of the 
velocity field the theory is superrenormalizable in the upper critical dimension fTsU , 
i.e. the number of superficially divergent diagrams is finite. In fact, it is sufficient to 
subtract the e-pole of the diagram shown in figure [I] to render the field theory finite 
at every order of the perturbation theory. The required renormalization factor can be 
obtained from the response function G\ } q = Go,i at one-loop order. A straightforward 
calculation yields 



J d d r exp(-ig,|r||)G 1 , (r,t) = exp(-^) 1 - - ^ _ (/2 f /2 + 0(Af 



e(l + e/2)(47r) 1 - 



(27) 



To subtract the e-pole in this function we introduce the renormalized diffusion constant 
Kr = Z^k, where 

2u 

Z = 1 - — and AJ = k r u^. (28) 

Here u is a renormalized coupling coefficient, A e = 1/ ((47r) 1 ~ e//2 (l + e/2)) is a geometrical 
factor and \x denotes an external momentum scale. 

Since the bare Green function is independent of fi its derivative at fixed bare 
parameters vanishes, i.e. 



_d_ 



G M ,Jv({r, *};*,/) = 0. (29) 



Figure 1. The only superficially divergent Feynman diagram of the field theory defined 
by the action S [equation (23)] with Ai = A2 = no = 0. The full line with an 
arrow represents the propagator (|2^) and the wavy line corresponds to the correlator 
fS(r± — r' ± ) of the velocity field. 



Expressing k and / in terms of Kr, u and fi we obtain for the renormalized Green 
function the renormalization group equation (RGE) 

\jidp - C{u)K R d Kn + P(u)d u ] G MtN ({r, t}; k r , u; = 0. (30) 

Due to superrenormalizability the Wilson functions 

H 

u = u ( K -e + 2u) (31) 



d/i 



In Z = 2u P(u) = U — 
o d ^ 





are exact at every order of the perturbation theory. Solving the RGE by characteristics 
one finds 

G M ,N({r, t}; k r , u; fx) = G MjN ({v, t}; Y(1)k r , u(l); /il) (32) 

where 

iilny(J) = -C(«(0) i^«(0 = W0) Y(l) = l u(l) = u. (33) 

For I — > the scale dependent coupling constant u(l) tends to the fixed point u± = e/2 
and y(Z) ~ XtZ _<E (with a non-universal scaling factor The asymptotic scaling 

behaviour of the Green functions follows from equation (|32| ) and dimensional analysis. 
The result reads (for I — ► 0) 

G M ,iv({r,t};^, M ;/i) ~ (/(F^)- 1 / 2 /^/ 2 )^ 

xG M ,v ({^(n^)- 1/2 / 1+£/2 r|| )/U /r ± ,/i 2 / 2 t} ; l) . (34) 

This equation shows that anomalous scaling dimensions of the variables and fields are 
given by 

r^-r 1 n ~r( 1+ ^ t~r 2 ^~i d+ *' 2 ~H~i d+e ' 2 - (35) 

The asymptotic form of the Green functions (34) depends on the arbitrary 
momentum scale fi and the non-universal amplitude Y*. A combination of these 
parameters occuring in (34) can be expressed by a quantity that is accessible to 
experiments or simulations |nj. Consider the random walk of a particle starting at 
time t = at the point r = 0. The mean square displacement of the particle in 
rn-direction averaged over the realizations of the velocity field reads 

X(t) 2 = J d rf rG lj0 (r,*)r 2 . (36) 

Using the one-loop result (27) and equation (34) with I = (/i 2 t) -1 / 2 it is easy to show 
that (for t — > oo) 

X(t) ~ at 1/2+e/4 with a = v^/^y*) 172 // 72 [l + 0(e 2 )] . (37) 

In the following sections we will always use the parameter a instead of Y*. Due to the 
strong anisotropy of the velocity field it is also possible to calculate X(t) directly from 
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equation (27) without application of the renormalization group. Since the n th order 
contribution to Gx$ is proportional to q 2n and X(t) 2 is proportional to the derivative of 
Gi t o with respect to q?< at q = the one-loop result is sufficient to obtain the function 
X{t) exactly: 

A '" )2 = M+ e(2 + f )(i)-<^ ' +,/2 ' (38) 
For e = 1 this result was first given by Bouchaud et al [[Tj|. 

For d = 3 (e = 0) an ultraviolet cut-off A is required to regularize the singularity in 
equation (|38|). In this case the mean sqare displacement is asymptotically given by 

X{tf ~ ^-tln(A 2 t). (39) 

27T 

4. Effective action below three dimensions 

The effect of the A2-vertex in equation (23) on the long-time behaviour of a reaction 
diffusion system without velocity field can be described by interactions located at the 
'time surface' t = 0. In reference |4]] it was shown by dimensional analysis that an 
interaction of the type 

d<V^-A m ,„ ^ (40) 
mini t=o 

is relevant below d m ^ n = 2(m + n) /(m + n — l) dimensions. (A^o and Ao,i are relevant in 
any dimension.) The critical dimension of A mn is changed if the motion of the particles 
is not purely diffusive, e.g. if the particles are subject to a linear shear flow ||11|| . In 
order to see how the values of d m ^ n are changed by the random velocity field consider the 
expansion of (4>{t)) in powers of Ai, A2 and /. For d < dj = 3 (d > df) the contributions 
to this expansion at a given order in Ai and A2 become more important (less important) 
for large t if we increase the order in /. Therefore the long-time behaviour of (4>{t)) for 
d > 3 is dominated by the diffusive motion of the particles and we may set / = 0. Below 
three dimensions we have to retain all orders in / to determine the asymptotic decay of 
the density. In this case the many-particle response functions ( p4"D play the same role 
as the Gaussian propagotor ( p5[) in the case / = 0. It is clear that it is not possible to 
sum the the power series in / exactly but the renormalization group allows us a partial 
summation and gives the correct scaling behaviour. 

We first determine the scaling dimension of A2 for a non-vanishing velocity field 
below three dimensions. Using the anomalous scaling dimensions given in (|35|) one finds 
A2 ~ i"2-d-t/2 _ j(i-d)/2^ - i e ^ j g j rre i eV ant for d > 1. In the same way we obtain 

A^~/ 2 A!0~/ 2 A^~/^ 2+e/2 A^0 ~ l d - 2+e/2 (41) 
and the dimension of the initial term 

\™ +n /\ mn ~ l0—m-n)(d+e/2)+2(m+n) _ ^2) 
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In d = 3 the vertices with m + n < 3 are relevant. For m + n > 3 the initial coupling 
A min is relevant below d m ^ n = (m + n + 3)/(m + n — 1) dimensions. However, we will 
show below that the only initial couplings which are generated by the coarse graining 
procedure employed in reference |3| are A 2>0 and A 0i2 (in addition to A 0) i = \/2n ), and 
these couplings are given by A 2 = — A 2 = n at every order in A 2 . 

Since the coupling coefficient A 2 is irrelevant for d > 1 one has to introduce a cut- 
off wave number A to avoid ultraviolet-divergencies in the perturbation series. The 
cut-off is required to regularize the diagrams shown in figure |2| which contribute to an 
effective coupling constant A 2je fr- The Green functions Gm,n considered in the previous 
section occur in these diagrams as subintegrals (represented by rectangles). As in the 
previous section we calculate these subintegrals without cut-off since they are ultraviolet- 
convergent for d < 3. After summation of the subdiagrams (using the renormalization 
group) the cut-off is introduced to perform the remaining integrations. If the cut-off 
wave number is not too large [an estimate based on equation (|38|) with t ~ A~ 2 gives 
A e <C f /(ne)} the scaling behaviour (34) still holds on the momentum scale of A. In this 
case the effective coupling has the scaling form A 2ie ff = A 2 Fa(A 2 A^ _1 ^ 2 /ct) . In the same 
way coarse graining leads to an effective coupling coefficient Ai e ff = AiF A (A 2 A( d_1 )/ 2 /cr) 
(with the same scaling function Fx). Evaluating the diagrams depicted in figure (|2|) 
by renormalization group improved pertubation theory it is possible to compute the 
coefficients of the Taylor expansion of F\ in an e-expansion. [Due to the scaling 
behaviour (34) of the response functions the diagrams are (infrared-) finite for d > 1 
in the limit of zero external momentum and frequency] However, this is not very 
interesting since the result strongly depends on the type of the cut-off. 

The initial coupling A 2 is more interesting since it is known to determine (at 
least for systems without random velocity field) the amplitude of (4>(t)) [§J. For the 



Figure 2. Contributions to the effective coupling coefficient A2, c ff to third order in A2. 
The hatched rectangles represent the Green functions Gm,n [equation where M 

(N) is the number of broken lines (full lines) running into the rectangle. The direction 
of each line is indicated by an arrow and reflects causality, i.e. each line points into 
the direction of larger time arguments. 
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calculation of A 2i0 it is convenient to shift the field (f> by the mean field density 

$ mf (t) = $o (1 + A^r 1 (43) 

with the initial value $o — v2no- The shifted action S[ip,ip; <j>, $(£) + 0] consists of a 
Gaussian part 

S G $, 4>, <f>; $(*)] = J & J d d r [$(dtip - A^ - Kdftfj) (44) 

+4> (d t <f> - A ±( p - Kd\<p + 2Ai$ mf (t)0)" 
and higher order terms summarized in 

S in S, i>; 4>, <P; = JdtJ d d r [Ai0(0 2 - ^ 2 ) + A 2 (2$ mf {t)<f> + $ mf (t) 2 )(0 2 - ft) (45) 

+A 2 (0 2 - ^ 2 )(0 2 - ft)] - f -J d d - l y (/ dt J dxtfd^ + 4>d^ 2 . 

Since Sq is independent of / the Gaussian (tree) propagators are the same as in 
reference ||, i.e. 

GV(q; t, t') = / d'Ve"^ ^(r, f$(0, f)>G = G(q; t - (46) 

G,(q; t, = | d rf re-^(0(r, t)0(O, f )>g = ( 1+^* ) G ^ t ~ ( 47 ) 

where G(q; i) is the propagator defined in equation (|26|). Corrections to the mean field 
density can now be calculated in a diagramatic expansion around the Gaussian action 
Sq treating the interactions in S int as perturbations. 

Following the lines taken in reference Q we replace the irrelevant parameter A 2 by 
effective initial couplings of the form fl4~0"|) . The diagrams contributing to A m n have 
(m + n) external legs each of which is associated with one of the vertices proportional 
to A 2 . Figure |3| shows the general form of the contributions to A 2 . Since the function 
®m{(t) 2 is damped on time scales which are large compared to (Ai$ ) 1 hs effect on 
the long-time behaviour of the field theory can be described by a 5-function. We may 
therefore write 



A 2 $mf(t) 2 ^ 2 



dt'$ mf (t 



>\2 



S (t)ft = ^ 5(t)ft (4? 



and obtain A 2j0 = 2A 2 $o/Ai + 0(A 2 ). In the same way higher order contributions to 
A 2j0 can be calculated by integrating the diagrams (b-e) in figure |3| with respect to the 
time argument associated with the leftmost vertex. However, the diagrams (b) and (c) 
cancel, and the same is true for the diagrams (d) and (e). In order to see this we write 
the contribution of the graph 0(b) in the form 

graph |(b) = / dt(-X 2 )f(t) (49) 
J 
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where t is the time argument carried by the leftmost vertex, and the function f(t) 
represents the subdiagrams in the hatched bubble and the loop integration. The same 



Therefore the contributions (b) and (c) cancel. Analogously one can show that the sum 
of the diagrams (d) and (e) vanishes. Since the same line of arguments applies if the 
hatched bubbles in figure |]have additional external legsf, all effective interactions A m>n 
with m + n > 2 are zero, and the only non-vanishing initial couplings are given by 
A2,o = —^-0,2 = 2X2^0/^1 — n o- In the appendix of reference Q a non-zero correction 
to A 2 ,o of the order O(\2 n o ) was obtained because the diagrams of the form (c) and (e) 
were incorrectly thought to be accounted for by taking Ai to A e fj. 

f Due to causality each diagram which contributes to the initial couplings contains at least one pair of 
i/j- or (j>- legs coming from a single A2 vertex. This is the vertex explicitcly shown in figure pi 



Figure 3. General form of the contributions to A2,o- The full and broken bold lines 
represent the propagators G^ and G^, respectively, and the zigzag line corresponds 
to the mean field density < & m f(t). At lowest order in A2 the bubbles in the graphs (b) 
and (c) contain only the graph (a) as a subdiagram. The vertex proportional to $ m f (t) 
in the diagrams (c) and (e) is generated by the shift of the field 4> [see equation (45)]. 



function f(t) occurs in the contribution [|(c) which is given by 



graph |(c) = / dm 2 $ mf (t) / dt'G^q = 0; t, *') Ai/(0 




(50) 




(a) 



(b) 



(c) 




(d) 



(e) 
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The cancellation of diagrams suggests that there should be a simpler way to calculate 
the initial coupling coefficients. In fact, the effective action can be derived without using 
Feynman diagrams. This formal derivation exploits the equation of motion 

D$, V; 0, fl-jr—expi-SW, V; 0, <t>\) = (51) 

and the corresponding equation with a functional derivative with respect to 0. The 
explicit form of the equations of motion (which are valid after insertion into averages) 
is given by 

d t tfj - A ± tfj - ncffa - 2A 2 ^(0 2 - V 2 ) - /(<9||V) / dt' J dx' (ipd^ + ^<f>) = (52) 
and 

d t <p - A ± - Kd\<p + A : (0 2 - ^ 2 ) + 2A 2 0(0 2 - ij 2 ) - %5(t) 

-f(d ]{ <f>) J dt' J dx' (i>d\^ + 0«9||0) = 0. (53) 

Averaging (53) one finds 

dMt) + A!(0(r,t) 2 -^(r,t) 2 ) = $ 5(t) (54) 

where $(t) = (4>(r, t)) is the exact density. (The equal time averages {(j)ip 2 ) and (4>4> 2 ) 
are zero due to the prepoint time discretization used in the derivation of the action.) 



Integration of equation fl54|) over t yields 



rod 

\J dt(0(r,t) 2 -^(r,t) 2 ) = $ o . (55) 
Jo 

The initial fluctuations generated by the irrelevant A 2 coupling can now be obtained by 
the replacement 

2A 2 (0 2 - V 2 ) — -> 2X 2 5(t) / dt'(0(r, t') 2 - ^(r, t') 2 ) = n 5(t) (56) 



o 

in the equations of motion ([52]) and (53). The effective equations of motion derived in 
this way are equivalent to the effective action with the initial vertices A m n calculated 
above. 

5. Asymptotic decay of the density 

In order to calculate the long-time behaviour of (<f){t)) we start from the effective action 
S e s = Srcac + S- mi , where 



+Ai0(0 2 - 4> 2 )} - ^ | d d ~ l y (j dt I dx(^|^ + 09||0)) 2 (57) 
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and 



$00 + -n (i) 2 - 2 ) 



t=Q 



(5? 



At this point it is convenient to reintroduce the Gaussian velocity field v(y) (with 
[^(y)] = and [v(y)v(y')] = fS(y — y') as in section ||) to replace S reac in favour of a 
f-dependent action which is local with repect to time. This new action is equivalent to 
the equations 



d t ip — A±ip — K,d 2 ijj + v(y)d\\ip = 

d t (j> - A ±( f) - KdR<f> + v(y)d»<f> + A(0 2 - ip 2 ) = 



(59) 
(60) 



where ip and are now real fields. The solutions ip, 4> have to be averaged with respect 
to both the initial conditions and the realizations of the velocity field. Since the initial 
state defined by Si n i and the distribution of v(y) are homogeneous the averages of the 
spatial derivatives in (59,60) vanish. Averaging equation (60) we find 



$(t)+A ($(t) 2 -(^r,t) 2 ) + C(t) 







where 



C(t) = ((0(r,t)-$(t)) 2 ). (61) 

To obtain the asymptotic solution of this equation one first has calculate (ip(r,t) 2 ) for 
large t. Later in this section we will also need higher moments {ip 2n ) with n = 1,2.... 
Since in equation (59) the dynamics of ip is decoupled from the reaction process the 
moments of ip may be calculated with the Green functions fl2~4|). In this way we get 

' -//,'! / <["/•', . . . / <["/•:. ; < '(r. /)-"(.-(r,.())- ...i -(r'... Or ) 

A=0 



(V(r,t)*> 



2 n n\ 

i const x (V^r^m 



(62) 



a 



where the scaling form (34) has been used. [A simple power counting as in reference [ H? 
shows that insertions of composite fields such as i\) 2n in equation (62) require no 
additional renormalizations.] 

To obtain the amplitude of the second moment (ip 2 ) at first order in e one has to 
calculate the diagrams shown in figure f|. The result reads in terms of the unrenormalized 
coupling constants 



(^(r,t) 2 ) = n «r 1/2 (87rf)- d / 2 



f^ 2 



(4vr; 



^-(i-3 ln2 + 0(e) ) +0(/2) 



Upon application of the renormalization group this becomes (for large t) 

V2 



(^(r,t) ) ~ n - 



-t 



-(d+3)/4 



l + -(31n2-l) + 0(e 2 ) 



(63) 



(64) 



We are now in a position to determine the asymptotic decay of $(£) in a similar way as 
in the case v(y) e0 1: 
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We first solve equation ([H]) for C{t) = 0, i.e. we are looking for a function <3? u (£) 
satisfying 

$ u (t) + A($ u (t) 2 -(^(r,t) 2 )) =0. (65) 

Since (ip 2 ) = est x t~( d + 3 )/ 4 for t — ► oo (with (d + 3)/4 < 2) the asymptotic solution of 
this equation is given by $ u (£) — Vest x t~( d + 3 )/ 8 . Using the positivity of C(t) one can 
show H that $ u (t) is an upper bound for $(t). To derive a lower bound $i(t) < $(£) 
we use the fact that a = (4> + ip)/ y/2 and b = (4> — ip)/ a/2 are real densities which satisfy 
the equations of motion 

d t a — Aj_a — nd 2 a + v(y)d\\a + V2Xab =0 

" " (66) 

d t b - A ± b - Kd\b + u(y)9||6 + V^Aafe = 0. 

For non-negative initial values the densities remain non-negative for all t. This implies 
2 - ip 2 = 2ab > and leads to the lower bound $i(t) = (|^|). 

Although the initial distribution of ip is Gaussian, the velocity field generates for 
t > higher cumulants in the distribution Vt of ip(r,t). We perform a cumulant 
expansion for this distribution to calculate (|V>|)- Since higher cumulants are generated 
by the disorder vertex / (or its renormalized counterpart v with fixed point value 
v* = e/2) this amounts to an e-expansion for Vt{ip). The Fourier transform Vt of 
the distribution with respect to ip can be written as 

where C2 n {t) denotes the connected part (cumulant) of the (2n) th moment of if), e.g. 



V t (h) = fjmW^ = exp (± (-l)^ 2n ) (67) 



c 2 (t) = (<Kr,t) 2 > c 4 (t) = (^(r,t) 4 ) -3c 2 (£) 2 etc. (68) 

For a Gaussian distribution all c m (t) with m > 2 vanish. Equation fl6"7|) can be used to 
calculate (\ip\) in an expansion in c 2n (t), n = 4, 6, . . .: 

(|^(r,t)|)=/ dm ^V t (h)e- ih * 



2tt 



Figure 4. Contributions to (?/> 2 ) at first order in e. The full and empty circles 
represent the initial vertex A2 j o and the field ip 2 , respectively. A solid line with an 
arrow corresponds to the Gaussian propagator (|2^). 
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'2c 2 (t) 



7T 



1 c 4 (t) 



+ 0(c 4 (t) 2 ,c 6 (t),...) 



4!c 2 (t) 2 ' w >-ow>--vj • (69) 

Due to the simple scaling C2 n (t) ~ £-«(d+3)/4 £ cum ulants for large £ the expression in 
the square brackets tends to a constant. The first correction to the Gaussian distribution 
comes from the four-point cumulant c±(t). The contributions to c 4 (i) at second order in / 
are shown in figure |5] (b-d). The first order (figure [| (a)) vanishes due to the x-derivative 
associated with the vertex /. This means that upon application of the renormalization 
group the amplitude of c 4 (t) at the fixed point is of the order e 2 . Analogously any 
cumulant C2 n {t) with n > 2 is of the order e n since a non- vanishing contribution to 
this function requires at least n vertices proportional to /. Therefore we only need the 
Gaussian part of the distribution to calculate (\if>\) at first order in e. The result is 

' V2n \ 1/2 



Ht) = 

Using the inequalities 



7T 87T 



) d / 2 a 



-(rf+3)/8 



l + -(31n2-l)+0(e 2 ) 



<({<P-W)<(<P 2 )-^ 2 ) 



(70) 



(71) 



in conjunction with the upper bound 3> u (t) as in one can show that the lower 
bound ( [TOD gives the exact long-time behaviour of the density, i.e. $(£) ~ $i(t) for 
t — > oo. 

In three dimensions the random velocity field gives a logarithmic contribution to 
the density $(t). Using equation one obtains 



1/2 



2tt 



-t 



-3/4 



(2/ln(A 2 t) 



-1/4 



(72) 



6. Summary and Discussion 



In this paper the effects of a quenched random shear flow on the long-time behaviour of 
the A+B-annihilation reaction have been studied. It is well-known that in dimensions 



Figure 5. Non-zero contributions to Ci(t) at second order in / (b-d). The first order 
diagram (a) vanishes. 




(a) (b) (c) (d) 
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d < 3 a random velocity field gives rise to enhanced diffusion and should, therefore, 
accelerate the reaction process. Similar to the case of purely diffusive particle motion 
the system segregates after a short time into regions of purely A or B particles. After 
this initial stage the density decay is governed by an effective Gaussian distribution for 
the initial density fluctuations. While a simple model based on anisotropic Levy walks 
already yields the correct exponent for the density decay one has to take many-particle 
correlations generated by the quenched randomness into account in order to calculate 
the amplitude. We have shown how this is possible in the framework of a systematic 
expansion in e = 3 — d and computed the first order term in this expansion. 

In the present paper it was assumed that the diffusion constants Da, Db of A and 
B particles are equal and that the hopping rate of both species depends in the same way 
on v (y). The more general case Da ^ Db with three different disorder couplings J'a, 
/b, Jab (instead of /) can be treated in a similar way at least as long as both species 
are mobile. For unequal diffusion constants the fields <fi and ip are coupled already in 
the Gaussian part of the action. Using similar arguments as in reference Q it can be 
shown that this effect changes only the amplitude of $(£). 

It would be interesting to check the results of this work by simulations. In order 
to observe the asymptotic long-time behaviour of a disordered system one has to 
perform configurational averages unless the system size is very large. In the case of 
a random shear flow the sample-to-sample fluctuations of the density depend on the 
linear size L y in the direction perpendicular to the velocity field. Since a diffusing 
particle covers a distance lj}(t) ~ {Dt) 1 / 2 during the time t the density measured at 
time t represents approximatively N ~ L y /ljj(t) independent configurations of the flow. 
Therefore, lack of self-averaging in a finite sample leads to a relative error of the order 

a pure system (without shear flow) the asymptotic power law 
can be observed at times of the order Dt ~ 10 5 Q. The above estimate shows that 
a random system with L y > 1.2 • 10 5 (or an equivalent number of smaller systems) is 
required to reduce the sample-to-sample fluctuations for Dt ~ 10 5 to less than 5%. 
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